
Available online at www.sciencedirect.com 

SCIE 



N c E ^o 1RE c T . PHYSICA 



ELSEVIER Physica E 29 (2005) 154-174 



www.elsevier.com/locate/physe 



Long and short time quantum dynamics: 
I. Between Green's functions and transport equations 

Vaclav Spicka a '*, Bedfich Velicky a ' b , Andela Kalvova a 

^Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Praha 8, Czech Republic 
^Faculty of Mathematics and Physics, Charles University, Ke Karlovu 2, Praha 2, Czech Republic. 

Available online 18 July 2005 



Abstract 

This paper summarizes the present views on construction of the electron quantum transport equations based on the 
Non-Equilibrium Green's Functions approach. The basic tool is the so-called Ansatz decoupling; one Ansatz family 
stems from the original Kadanoff-Baym Ansatz and is suitable for extending the quasi-particle picture of the Landau 
theory out of equilibrium. The other family based on the Generalized KB Ansatz is appropriate for short time 
transients. The physical and formal context of the Ansatzes is analyzed; the most important question explored is the 
status of the Reconstruction Theorems, reducing a full description of a non-equilibrium system to a dynamic theory in 
terms of one-particle quantities. A comparison is made between the Time Dependent Density Functional Theory and 
the properly renormalized NGF formalism. There is a close relationship between both formalisms. The Reconstruction 
Theorems form a general basis for obtaining improved quantum transport equations. 
© 2005 Elsevier B.V. All rights reserved. 
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1. Introduction 

This paper is the first one of three contributions 
to this volume ([1], hereafter called Paper II, [2], 
hereafter called Paper III) devoted to quantum 
transport out of equilibrium, that is to say, beyond 
the linear response. We concentrate on a specific 
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question of generating reliable transport equations 
for electrons in extended systems driven by external 
fields under various conditions which might be 
broadly divided into "slow" and "fast" phenomena, 
represented by quasi-stationary quantum transport, 
and optical or transient processes, respectively. We 
need to describe both regimes properly. As the 
underlying general description of non-equilibrium 
systems, we will use the Non-equilibrium Green's 
Functions (NGF) [3-9] - double-time correlation 
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functions which include full quantum sub-dynamics 
of interactions in the presence of arbitrary external 
fields and/or in highly non-equilibrium transient 
states. The NGF can be used to describe actual 
systems directly, but we presently discuss the 
methods of an approximate reduction of the theory 
to equations for single-time quantities, quantum 
distribution functions. The approximations are 
based on hierarchy of characteristic times of a given 
system and they all have "Ansatz" in their name 
with reference to the original, so-called Kadanoff- 
Baym Ansatz (KBA) [3,8]. 

In agreement with the above division of 
phenomena to slow and fast, the basic formalism 
forks into two specific lines of developing single- 
time evolution equations. The common "intuitive" 
tool to describe transport experiments (long time 
behavior) is a quantum generalization of the 
Boltzmann equation (BE) [3-12]. It is based on 
the concept of a quasi-particle distribution func- 
tion and on the quasi-classical approximation 
suitable if everything varies smoothly in space 
and time. In the NGF context, the mild time 
dependence, in particular, is a necessary condition 
for the validity of the KBA, which has been 
devised to obtain the quantum generalizations of 
the Boltzmann equation. We will call this type of 
evolution equations the Quantum Kinetic Equa- 
tions (QKE). Further progress in this direction led 
to an extension of the Ansatz incorporating off- 
shell propagation in an effective manner. The 
resulting modified scheme has been used to obtain 
the quantum BE with the corresponding correc- 
tions [6,13-15]. 

If a smooth time dependence cannot be 
assumed, a more complete quantum treatment 
becomes necessary to describe the short time 
(transient) and high-frequency processes. For 
states far from equilibrium, a theory interpolating 
adequately between the fully quantum linear 
response (...the standard Kubo formula, whose 
use is limited to weak external fields acting on an 
equilibrium system) and the quasi-classical high- 
field BE is needed. Such a theory may be based on 
the quantum Generalized Master Equations 
(GME) [8,9,16,12] with a strictly causal structure. 
These equations abandon the quasi-particle dis- 
tribution concept and describe the evolution of the 



density matrix (or Wigner function) for true 
particles. The tool for generating the GME within 
the NGF environment is a correspondingly 
structured Ansatz usually called Generalized KB 
Ansatz (GKBA) [17-20,8,9]. Also this Ansatz 
can be further developed and refined to obtain 
an improved GME with a wider range of applic- 
ability [8,9]. 

It is remarkable that more than 40 years since 
the inception of the KBA, and almost 20 years 
since its modification to the GKBA, this area of 
non-equilibrium physics is still active and devel- 
oping, as witnessed by the number of papers and 
of meetings devoted to it [21-23]. The questions 
addressed range from the practice of employing 
the quantum transport equations to the justifica- 
tion, extension and refinement of the Ansatz-based 
bridge between the NGF and other quantum 
transport approaches. These problems are the 
subject matter of the present paper series. 

It should be stressed that similar questions and 
tasks emerge in very small metal and semiconduc- 
tor systems. Quantum coherence plays an essential 
role in the behavior of such systems. A natural 
question emerges to how the behavior of such 
systems corresponds to the "classical" way of 
description known from classical thermodynamics 
and statistical physics. We attempt to describe a 
unifying approach to non-equilibrium dynamics, 
which, while based on studies of extended systems, 
is also well suited for systems with a reduced 
dimensionality. 

Summary of Papers I, II, III 

This Paper I and the adjoining Papers II, III are 
not intended as an exhaustive review. They rather 
try to stress the basic issues from our point of view. 
This is especially true of this first one. The subject 
matter is divided roughly as follows. In this paper, 
we present an overview of the basic concepts and 
of the formal means used. The paper culminates in 
the last section devoted to the reconstruction 
problem: to what extent can the non-equilibrium 
dynamics of the many-particle system be described 
in terms of the evolution of the one-particle 
distribution? Paper II is devoted to the case of 
the quasi-particle transport in quasi-classical 
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environment, as described by the QKE. The 
Kadanoff-Baym Ansatz and its extensions are 
analyzed in full detail. Finally, Paper III studies 
the GME approach to strongly non-equilibrium 
processes. This is connected with the physics, 
improvements and generalizations of the GKBA. 
Special attention is paid to the proper treatment of 
initial conditions, which are crucial for the early 
evolution of transients. 

2. Definitions and symbols 

The average values, for all additive (single- 
particle) observables A(x,p), are given by the 
single particle reduced density matrix: 

AUt) = tr(p(t)A) , (1) 
p(x,y, t) = TR(p^(y, W(x, t)), 

EE <<lAW)<K*,0>> • ( } 

Here p is the full many-body statistical operator 
(typically, a grand canonical density matrix) and 
i/^ and xj/ are field operators for particles (electrons 
throughout this article), all in the Heisenberg 
picture. 

In particular, the local particle density n of 
electrons is given by 

n(x, t) = p(x, x, t) . (3) 

For a single band e(k), the local current density j 
equals 

^•' )= ^(^ )+ S ( - , ^ ) )^ ,) U- 

(4) 

A good starting point for the real-time NGF 
approach is to extend, following Kadanoff and 
Baym [3], the single- time reduced electron density 
matrix p to a double-time function, the so-called 
electron correlation function 

g < (\,2)= « ^(2)^(1) », (5) 

where, instead of x, t we introduce cumulative 
variables denoted by numbers, l=x\,t\, a = 

The single-electron reduced density matrix is 
then the time-diagonal element of this correlation 



function 

p(x u x 2 ,t) = g < (h2)\ h2=r (6) 

A complete set of correlation functions is 
obtained by introducing a complementary func- 
tion, the correlation function for holes 

^(1,2)= « ^(1)^(2) ». (7) 

For these two functions g < , g > , Kadanoff and 
Baym developed a closed set of coupled dynamic 
equations, the KB formalism [3,8]. In parallel, 
Keldysh created a basically equivalent matrix 
formulation of the NGF approach [24,25,5]. See 
Section 7. 

Remark on Conventions: This section intro- 
duces most of the writing conventions used in this 
paper. There is no special symbol for single- 
particle operators, while many-body quantities 
have a cap. Single-particle trace ...tr, many-particle 
trace ...TR. Commutator ...[...]_, anti-commutator 
...[...] + .We put h = 1, so that there is no distinction 
between p = hk and k. The dispersion law e{k) = 
£k need not be quadratic, a general band shape is 
admitted, so that the kinetic energy of the 
electrons is T = ^ \k)eu{k\. 

3. Physical concepts 

While a true historical overview would be 
beyond the scope of this paper, it will be impor- 
tant to recall several fundamental ideas which 
run through the development of the transport 
theory. As will be seen, the physical principles 
of the contemporary NGF/Ansatz approach are 
not different from those which were on the minds 
of the founders of non-equilibrium statistical 
physics. The issues remain, the understanding 
and attitudes here changed. We choose just three 
topics. 

I. The central idea of the transport theory has 
always been a reduced description of the system. 
The relevant information about a gas of particles 
was contained in the one-particle distribution 
function f(x,p, t). The essential step was to search 
for a closed equation governing this function, that 
is for the transport equation. This was only 
possible, if, say, the description of binary collisions 
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in the gas was given in terms of the function /. As 
a systematic program, this seems to appear first in 
the work of Chapman and Enskog [26,27], but the 
principle dates back to Boltzmann himself [28-32], 
with his Stosszahlansatz for the binary distribution 
f n expressed as fn^fi x f 2- This type of factor- 
ization, or, more generally, decoupling of higher 
correlation functions was systematically developed 
in the BBGKY technique. We have to say more 
about Bogolyubov, to whom we owe the postulate 
that in a chaotized many -body system all higher 
particle correlation functions become a functional of 
the distribution function [33-35,9,16,12]. Thus 



relaxation of spatial inhomogeneities. In "normal" 
situations, the three times obey the inequality 



chaot. 

fn — ► WxJi] ■ 



(8) 



The functional has yet to be specified, but this is of 
secondary importance. As will be described later, 
the KB Ansatz also singles out the one-particle 
distribution as the determining characteristic of 
the non-equilibrium system. This brings us to the 
reconstruction problem: under which conditions the 
full description of the many-body interacting 
system can be built up from the knowledge of 
single-particle characteristics? This seemingly out- 
rageous question was seriously treated in several 
contexts. In Section 8 we will discuss two of them: 
the time-dependent density functional theory 
(TDDFT), and the generalized Kadanoff-Baym 
Ansatz related techniques. 

II. The second crucial notion is the hierarchy of 
characteristics times. There are three intrinsic times 
related to a many-particle interacting system. In a 
reminiscence of a non-dilute gas, they are often 
identified as the collision duration time x c = a/v, 
the collision time x r = £/v and the hydrodynamic 
time Th = L/v. Here, v is the average thermal 
velocity for a classical gas, the Fermi velocity 
for a degenerate Fermi gas. The characteristic 
lengths are: a the interaction potential range 
(particle size), £ the mean free path (mean inter- 
particle distance), and L the characteristic length 
of spatial inhomogeneities in the system. In the 
modern interpretation, the three times are: x c • • • 
the chaotization time characterizing the decay of 
correlations, x r • • • the relaxation time characteriz- 
ing the thermalization of the system (local relaxa- 
tion) and, finally, Xh • • • characterizes the process of 



T c < T r < X h 



(9) 



separating the chaotization stage, the kinetic stage, 
and the hydrodynamic stage. This structuring of 
the spontaneous return to equilibrium was also 
introduced by Bogolyubov; he postulated Eq. (8) 
only for times later than the chaotization time x c . 
The Boltzmann equation proper corresponds to 
the limit 







X c <^ T r <^ Th 00 



(10) 



We will be mostly concerned with the opposite, 
more realistic case, when the distinction between x c 
and x r will be less sharp. Further, the time range of 
interest will be specified by the external fields. For 
example, an optical pulse is characterized by its 
duration, the ground period of the signal and its 
Rabi period measuring the pulse strength. These 
times should be compared with the intrinsic times 
of the system. This will, finally, specify the 
situation and the necessary version of transport 
theory used. 

III. The last basic concept is that of the quasi- 
particles. With this, we, as a matter of fact, move 
over to the quantum realm. There are several 
streams merging into the generalized notion of a 
quasi-particle. Firstly, the polaron, an electron 
dragging along a cloud of lattice polarization. In 
the quantum field language, the electron is dressed 
by virtual phonons. This compound object has 
some features characteristic of a particle, like a 
dispersion law (renormalized by the self-energy), 
which has an operational meaning in experiments. 
The self-energy is typically complex, and this leads 
to the finite life-time x of the quasi-particles, 
closely related to the transport relaxation time x r . 
Secondly, and even more to the point, in a non- 
dilute system of interacting particles, their indivi- 
duality is suppressed by mutual correlations and 
the use of a transport equation seems to be useless. 
However, the weakly excited states may appear to 
mimic a gas of weakly interacting quasi-particles. 
This was at the bottom of the Landau theory of 
the Fermi liquid [36], whose part was a proper 
adaptation of the Boltzmann equation [25,37]. 
This approximate, but highly precise description of 
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the Fermion systems had, of course, a number of 
predecessors, like the Sommerfeld electrons in 
simple metals, and parallels, like the quasi- 
particles in nuclei, where the self-energy has been 
originally introduced under the name of the 
optical potential. It might seem that the quantum 
transport equations for strongly interacting sys- 
tems will all deal with the quasi-particles. This is, 
to some extent, true. However, the quasi-particles 
are vulnerable and elusive objects and cease to 
exist under some harsher conditions, like beyond 
the quasi-classical regime, for strong and/or 
transient disturbances, or if the system itself is 
not favoring their existence. This can sometimes be 
judged by the Landau-Peierls criterion [38]. Take 
the gas of quasi-particles with energies around the 
Fermi level and a lifetime t due to impurity 
scattering. The criterion reads: if t — h/E?, the 
BE-like transport theory is not applicable. Now 
h/E F can be interpreted as the "quasi-particle 
formation time" tq. Ifr ~ tq, the quasi-particle 
decays before having formed. On the other hand, % 
is closely related to r r , whereas tq appears to play 
the role of t c . The Landau-Peierls criterion thus 
states that if the quasi-particles do not form, then 
u r becomes comparable to t c and the condition 
Eq. (10) is not obeyed. 



4. Origins: Boltzmann equation for classical and 
quantum particles 

A prototype of all kinetic equations, the 
Boltzmann equation of classical physics 
[28-30,26,27,31,32,16,12] reads 



dt 



■drift [/■(*)] = Iin[f(t)]-Iout[f(t)]> 



(11) 



where the drift term is the classical Poisson bracket 

dsdf dsdf 



drift [/Xr)] = 



dk dr dr dk 



(12) 



This is, in principle, an equation for a single- 
particle distribution function f(k,r,t) in phase 
space, representing balance between the drift of 
particles (with energy s = + U(r 9 1) and velocity 
the external force field 



m 



dk 8 ~ dk Ck 

— S- s = — S- U), as given by the left-hand side of 



the BE, and the irreversible evolution due to 
collisions described by the scattering integrals I m 
and 7^ on the right-hand side. In these integrals, 
the collisions are approximated as instant rando- 
mizing events, so that the BE is Markovian. This is 
in agreement with the first inequality in Eq. (10). 
The other inequality demands all inhomogeneities 
in the system, including external fields, to be 
smooth enough to allow sufficient time for local 
equilibration. 

The BE has been used for many classical systems 
and it was extended very early also to the trans- 
port by quantized particles, in particular, by 
electrons in metals. The quantum effects were 
incorporated mainly in the scattering integrals 
where the collision rates were calculated by the 
Fermi Golden Rule and the exclusion principle 
was taken into account by means of the Pauli 
blocking factors. 

The average values of additive observables were 
calculated from formulas taken over from classical 
kinetic theory of gases. Thus, the particle density 
and the particle current density were 



n(r, t) 



dk 

(2n)- 



;f(k,r,t), 



... f dk de 
j(r, t) = / ^ — -f (k, r, t). 



(13) 
(14) 



Now, we depart from the BE in the direction 
towards quantum dynamic equations far from 
equilibrium. There are several extensions of the BE 
that we will first explore separately: 

(i) From particles to quasi-particles. The interac- 
tion energy will no longer be negligible, but 
will be renormalized in the quasi-particle 
transformation. 

(ii) From near-equilibrium to "arbitrary" non- 
equilibrium. This will be reflected in the 
memory of the system and taken into account 
by the non-Markovian GME. 



5. Quasi-particles and transport 

In dense systems, the particle interactions 
cannot be reduced to rare collisions randomizing 
the motion of otherwise free particles. Quantum 
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statistics for such systems at near-equilibrium are 
often well described by the Landau theory of 
quasi-particles. In this theory, weakly excited 
states of the interacting particles are described as 
a gas of quasi-particles with energy 

s(k, r, t) = e k + U Q ff(r, t) + a(s, k, r, t), (15) 

where £/ e ff m ay include the mean-field part of the 
interaction and the self-energy a describes the 
mass renormalization. The quasi-particles are 
coupled by a weak residual interaction. 

5.7. Quantum kinetic equation for quasi-particles 

These quasi-particle features have consequences 
for transport properties, which are described by 
the Landau kinetic equation. It has the exact 
structure of BE (11). However, the energy entering 
the drift term is now the quasi-particle energy s. 
Scattering integrals I m and 7 0Ut contain scattering 
rates calculated again by the Fermi Golden Rule, 
but with the residual interactions reflecting that 
they are reduced by the many-particle wave-func- 
tion renormalization. Finally, the function f(k, r, t) 
is the quasi-particle distribution function. 

Because of its intuitive character, the "Boltz- 
mann" equation for quasi-particles provided many 
insights into the behavior of non-equilibrium 
many-particle systems. In particular, it offered a 
retrospective explanation of the over-successful 
Sommerfeld model of metals with non-interacting 
electrons responding to external fields. While 
various modifications of the BE-like approach 
have been applied successfully to many systems 
and situations, there are also many physical cases, 
in which it is bound to fail. The limitations of the 
Boltzmann equation include: (1) The BE well 
describes quantum systems if the quasi-particle 
picture is justified. This may no longer be valid in 
highly non-equilibrium quantum systems. (2) It is 
based on the intuitive idea of instant collisions 
between quasi-particles, which requires the colli- 
sion duration time to be very short. (3) At the 
same time, the BE will not be well suited to 
describe systems with abrupt changes in space, 
including small structures, where quantum effects 
are essential and the quasi-classical approach 
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hidden behind the BE is far from being sufficient 
[8,9,39,40]. 

5.2. Observables and quantum distribution function 

Relations between the quasi-particle distribu- 
tion function and the expectation values of 
observables are more complicated than (13) and 
(14). For example, the expression for current 
density must incorporate the back-flow accompa- 
nying the motion of quasi-particles. Rules for 
evaluating observables have been outset from an 
integral part of the Landau equilibrium theory. A 
natural question emerges: which is the relation 
between the quasi-particle distribution function / 
and the expectation values out of equilibrium? 

Our task is then to find a functional relation 

{f, A}^A Av = tr(pA), 

which is equivalent to finding a functional p[f] 
associating the single-particle density matrix with 
any given quasi-particle distribution. To relate the 
reduced density matrix p with the distribution 
function / (k, r, t) it is convenient to express p in the 
(k, r) representation. This can be done in different 
ways, but we choose the so-called Wigner repre- 
sentation 

p(k, r,t) = J dxQ~ lkx p (r + | , r - | , /) . (16) 

The reduced density p(k,r,t) in the Wigner 
representation is commonly known by the name 
of Wigner distribution. It can be viewed as a 
distribution of electrons in the phase space and the 
expectation values 

n(r,t)= [ -^p(k,r,t), (17) 

J (271) 

»■»- J m 

have the same form as in the Boltzmann theory. 

These suggestive properties of the Wigner 
distribution should not lead us to the conclusion 
that this function is the "right" quantum general- 
ization of the classical distribution function, and 
that the proper quantum generalization of the 
Boltzmann equation will be a kinetic equation 
of the Boltzmann form (11) for the Wigner 
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distribution function. We can see that this is not 
the case from the following argument given 
already by Landau. 

Consider a homogeneous system. In equilibrium, 
the BE-like kinetic equation for electrons is solved 
by the Fermi-Dirac function regardless of the 
interactions in the system. Thus, at zero tempera- 
ture, the distribution / jumps from 1 to at the 
Fermi level. In contrast, the Wigner distribution 
describes the occupation numbers of true particles 
and it differs from the Fermi-Dirac function by 
depletion of the momentum states below the Fermi 
level. These missing states emerge as states above 
the Fermi level, where they form the so-called high- 
momenta tails of the Wigner function. The step of 
the distribution p is accordingly reduced to 1 — 2z, 
where z is the renormalization constant for quasi- 
particles at the Fermi level. 

We may summarize that the quasi-particle 
kinetic theory consists of two steps. First, the 
quantum kinetic equation is solved for the quasi- 
particle distribution. Second, the true quantum 
particle distribution is constructed by means of the 
functional p[f] and the expectation values of the 
observables are calculated. 



6. Generalized Master Equations: beyond the 
quantum kinetic equations 

As already mentioned, quantum kinetic equa- 
tions of the Boltzmann type have a restricted range 
of validity and these limits can be transgressed 
only by resorting to a more general framework, 
permitting, at least in principle, to start from a 
fully quantum description and work directly 
with an equation for the reduced density matrix 
(or the equivalent Wigner function), the so 
called quantum GME. There are several 
approaches and approximations leading to the 
GME [41-44,12,8,9], which has a general form 

- drift[p(0] = interaction term, 

^ + i[ r+y eff (0 ,p(Q]_ = J^dtFipit)]. 

(19) 



This equation is shown in an entirely symbolic 
form on the first line, and it is not actually explicit 
on the second one either. We will discuss the 
detailed structure of this type of equation below. 
Here, we only focus on its most salient features. 
The GME is a closed equation for p. It is non- 
Markovian, because the interaction term is non- 
local in time. 

At the end of Section 5.1, we pointed out three 
important limiting factors for the use of the 
quantum kinetic equations. These factors will now 
be discussed from the point of view of the GME. 

6.1. GME is an equation for p 

Eq. (19) is a dynamical equation for the 
distribution of particles instead of quasi-particles 
as was the case of the BE. Therefore, it does not 
hinge on the use of quasi-particles, and it is free 
from the physical limitations necessary for intro- 
ducing the gas of quasi-particles. In particular, the 
system may evolve under conditions which do not 
permit the quasi-particle states to consolidate. This 
does not preclude using some of the quasi-particle 
features where applicable. Thus, the left-hand side 
of the equation describes bare particles drifting 
under the influence of the effective field. All other 
features involving interactions are included in the 
right-hand side. It may be possible to transfer 
parts of the r.h.s to the drift term and achieve 
evolution in re-normalized bands. 

6.2. Interaction term 

We prefer to call the right-hand side of (19) on 
interaction term rather than scattering integrals. In 
principle, the GME is exact, so that the interaction 
term must incorporate all of the sub-dynamics of p 
reflecting not only particle collisions, but also the 
short time dynamics, off-shell propagation and 
coherence between the collisions and with the 
external fields, multiparticle correlations, gradual 
saturation of the scattering rates after the onset of 
a non-equilibrium process, etc. 

It is remarkable that all this rich physics can be 
absorbed in an interaction term depending only on 
the one-particle density matrix. This is made 
possible by the subtle memory effect reflected in 
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the time integration over the full depth of the past 
in (19). Clearly, this form of the interaction term 
agrees with the Bogolyubov postulate quoted 
above in Section 3, at least for times beyond x c 
after the onset of the process. We will return to the 
related Reconstruction Theorems in Section 8. 

6.3. Quasi-classical expansion 

In the GME (19), neither the drift term, nor the 
interaction integral, are restricted to smooth 
variation of the fields and distributions in space. 
The drift term is given by the quantum Poisson 
bracket rather than by the classical one appearing 
in (12). For smooth functions, it can be quasiclas- 
sically expanded; the formal expansion parameter 
is h 2 . For the Wigner distribution p defined in 
Eq. (16), we obtain (writing, in the way of 
exception, the Planck constant explicitly): 



drift[p(0] = Pr 1 

3 



[T+U e{f (t),p(t)]_ 
dU ef{ dp 



dr dp 



3 3 e 8 3 p 
dp 3, 8r 3 



9 3 ^eff 9 3 p 

8r 3 dp 3 



+ . 



(20) 

This expansion, suited for comparison of the 
GME in the quasiclassical limit with the corre- 
sponding kinetic equation is, in fact, a Fourier 
transformed expansion around the space diagonal 

X\ — X2. 

7. NGF approach to non-equilibrium systems 

This section only serves to introduce the 
concepts and relations concerning the NGF as 
they will be needed for our purposes, namely to 
derive and analyze the non-equilibrium dynamic 
equations of either the BE or the GME type. The 
formalism is commonly known today and there 
are a number of texts suitable for Refs. [3-9]. 

7.1. Green's functions on a time contour 

In a short span of time, Schwinger and Martin 
[45,46,21], Kadanoff and Baym [3,22] and Keldysh 



[24,23] created the field theoretical perturbation 
scheme for the non-equilibrium Green's functions 
starting from causal Green's functions 



G(l, 2) = -i« T c [^(l)^(2)] » 



(21) 



with the time-ordering operator T c acting on a 
time contour C . Here, in Fig. 1 we illustrate the 
Schwinger contour, usually referred to as the 
Keldysh contour. 

The rightmost turning point, just like the 
terminal point is usually sent to + or — infinity; 
we prefer to keep to finite to make contact with the 
problem of the initial conditions otherwise elimi- 
nated by the limit to — oo. For given times t\, ti, 
there are four different orderings on the C 
contour, as shown in the figure, and an equivalent 
formalism working with a matrix Green function 
of real times running from — oo to +oo, 



D 



G c G< 
G > & 



(22) 



can be established. Here, G c = — i^ T[^i/r]> and 
G = — i<^ J[i/^]> are the usual causal and anti- 
causal GF's of real times, while 



(7>(l,2) = -i« iKl)^(2) » 
G < (l,2) = +i«^ ( 2) lA (i)» 



(23) 



differ only by factors of i from g K , g > as defined 
in Eqs. (5) and (7). Sometime, it may be convenient 
to embed the C contour in the complex time plane 
in order to use the procedure of analytical 
continuation to link the contour with the Matsu- 
bara imaginary line [3]. In general, however, the 
contour has an entirely symbolic meaning. 

Only two of the four functions (22) are 
independent. Kadanoff and Baym took for these 
two functions G > ,G < and developed their KB 
formalism [3]. In parallel, Keldysh created a 
basically equivalent matrix formulation of the 



,G > _ 



Fig. 1. NGF time contour. 
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NGF approach. There is quite a choice among the 
variants of these two approaches [3-9]. To retain 
the matrix structure of the equations, a slightly 
redundant set of three functions is needed. We will 
use a modified Keldysh formulation suggested by 
Langreth and Wilkins [47,48]. It is obtained by a 
linear transformation 



manipulations. As a basic and typical example 
consider the product [£G] C on the contour. In the 



M 



Gh+M<7 Z GM 1 = G 
with 

1 0" 

1 1 

Result: 

G R G< 

G A 




-1 



(24) 



(25) 



(26) 



The transformed set consists of Green's functions 
including G*^ and the retarded and advanced 
propagators G R , G A defined by the identities 

G C (1,2)-G < (1,2)=G R (1,2), 
G C (1,2)-G > (1,2)= G A (1,2), 
G C (1,2)-G < (1,2) = -G A (1,2), 
G C (1,2)-G > (1,2) = -G R (1,2). 



(27) 



Explicitly, 

G R (1,2) = 6{h - t 2 )(G > (1,2) - G < {\,2)), (28) 

G A (1,2)= -6{t 2 -h)(G > {\,2)-G < {\,2))- (29) 

Keldysh used the set G R , G< + G> , G A instead 
[24]. 

The fourth GF is expressed in terms of the 
triplet through the obvious but important "spec- 
tral identity" 

G R (1,2)- G A (1,2) = G > (1,2)- G < (\,2) 

= -L4(l,2). (30) 

The triplet is overcomplete, as there holds the 
conjugation G R (1, 2) = [G A (2, 1)]* (any of the 
three functions G R , G A , the spectral density, is 
fully equivalent to the other two). A further 
reduction to a single GF is only possible in 
equilibrium. 

The set G R , G < , G A arranged to the upper 
triangular matrix G is extremely useful in algebraic 



matrix form, it is 


given by 








"G R G<" 


IG = 










Z A 




G A . 



Z R G R 




Z R G < + 2' 
Z A G A 



(31) 



These are the "Langreth- Wilkins" rules in a 
compact form [47,48]: the retarded (1,1) compo- 
nent contains only R-factors, similar to the 
advanced component (2,2). The < ("less") (1,2) 
component is a sum of all possible products with 
factors in the causal order R •••<••• A. These 
rules hold for any number of factors (chain rule). 
It is important that the less component of the 
product is necessarily linear in the less quantities. 
In an explicit form, all products are integral over 
the internal variables, like 

[G R Z < ](1,2)= J d3G R (l,3)I < (3,2) 
= G R (l,3)i: < (3,2). 

7.2. Explicit construction of the NGF 

Up to now, the procedure of introducing the 
NGF has been formal, but wholly general. It offers 
no practical means of constructing an actual GF 
for a specific situation. We have to recall that the 
<l • • • y average means, in fact, TR(p • • •), where 
the many-body initial state p at t = to has so far 
been arbitrary. If, say, a perturbation expansion in 
the powers of the interaction is intended, the initial 
state should belong to an unperturbed Hamilto- 
nian quadratic in the quantum fields; this is barely 
the case for any finite to, and we have to resort to 
the usual limit to —> — oo associated with the 
assumption of the "decay of initial correlations", 
another of the Bogolyubov tenets. This is the 
Keldysh way [24,5,25]. Alternatively, the trajectory 
may be bent along the Matsubara interval {to, to — 
iP) in the complex plane, but again with to —> —oo. 
This is the Kadanoff-Baym variant [3]. All other, 
more general situations amount to dealing with the 
so-called problem of correlated initial conditions, 
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which has been treated by a number of authors 
[49-63]. Less attention has been paid so far to the 
connection between the initial conditions and the 
construction of the transport equations [55,56,60]. 
This is natural for a steady transport: the initial 
correlations die out while the steady regime sets 
on. The situation may dramatically differ for fast 
transients and most of Paper III will deal with the 
related questions. 

Here, we prefer to concentrate on the basic 
Ansatz formulation and on the related reconstruc- 
tion problem. Therefore, we will only deal with the 
physical situations, in which the initial conditions 
pose no additional problems. Specifically, we 
restrict our considerations to the class of the so- 
called switch-on states, that is, we start from an 
equilibrated state disturbed subsequently by an 
external field. 

In this equilibrium starting state including the 
interactions, there exists an unperturbed state, as 
required by the principle of weakening of initial 
correlations. The two states could be linked by an 
equilibrium perturbation scheme. It is preferable 
to first act on this unperturbed state by the 
(effective) external field. This permits to introduce 
the reference unperturbed field-dependent NGF 
matrix Go incorporating already the required 
boundary conditions. Finally, we turn to the 
interactions by a non-equilibrium technique. The 
resulting full NGF defines the self-energy matrix 
from the Dyson equation 

G = G + G 1G, (32) 

which can be handled by the LW rules, etc. 

To close the equations, it is necessary to specify 
the (approximate) self-energies I ,< , Z R,A . For a 
consistent transport theory, it is important to 
express these in terms of Green's functions 
themselves, that is, to use a self-consistent 
approximation. To be physically correct, the 
approximation should further obey the basic 
requirements, like the symmetry laws, conserva- 
tion laws, causal structure, etc. 

We will not discuss specific approximations for 
the self-energies, assuming that the self-energies in 
question are physically consistent, so that the 
related Green's functions can serve as a basis for 
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constructing transport equations for quantum 
distribution functions. 

7.3. Dyson equations for propagators; self-energies 

By the LW rules, the non-equilibrium propaga- 
tors G R,A satisfy the usual Dyson equations. 
Written in the symmetric differential form 

{G R,A } _1 = Go 1 - Z R < A , (33) 

they serve, in fact, as defining identities for the 
respective self-energies. We have to specify the free 
Green's functions G R,A . As indicated in (33), they 
have a common inverse Gq . In agreement with 
(32), it is defined as incorporating the effective 
external field augmented by the mean-field con- 
tributions to the potential energy, 

Go _1 a,2) = S(l - 2)(i J- - e(j V 2 ) - U^(x 2 ,t 2 )^. 

(34) 

The self-energy I is cleared of the mean field, that 
is of a "singular" time-local component, by this 
definition. There are altogether four combinations 
as differential or integral Dyson equations corre- 
sponding to (33); we quote one of them for 
reference: 

G " 1 G R = / + Z R G R , (35) 
G R = G R + G R Z R G R . (36) 

Here, 1(1, 2) = 3(1 — 2). The effective external 
fields and the boundary conditions, that is, the 
jump conditions at equal times (from zero to — i 
for G R , from +i to zero for G A ) enter the problem 
through the free propagators G R ' A . The integral 
Dyson equation serves the important purpose 
of transmiting these properties on the full pro- 
pagators. 

7.4. Dyson equation for G < 

Equation of motion for the correlation function 
G < will play a central role in this paper. It is 
convenient to start from its differential form. To 
this end, we left-multiply Eq. (32) by Gq 1 and then 
take the less component. Employing (33), we 
finally obtain the first of the following two 
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conjugate integro-differential equations of motion 
for G < , valid for times t\ > to, t 2 > to'. 

{G R }- 1 G < = Z < G A , (37) 

G < {G A }- 1 = G R I < . (38) 

These equations are fully renormalized. We have 
to treat the initial/boundary condition with great 
care. Either we perform the limit —> —oo, or we 
decide to keep to finite. Then, at variance with the 
propagators, the initial condition does not just 
reduce to a unit jump at equal times. Instead, the 
conditions at an initial time to have to be specified 
(do not mix this unfortunate "0" with that of the 
free propagators). In the simplest case, taken into 
consideration by Keldysh and corresponding 
neatly with the transport equation approach, the 
initial condition is assumed as uncorrelated and 
reduces to the knowledge of the one-electron 
density matrix 

P(*0) = Pinif ( 39 ) 

In that case, G < has the form, frequently called the 
Dyson-Keldysh Equation, 

G< = G R p< lt G A + G R I < G A , (40) 
where 

p 4(1, 2) = id{h - t )S(t 2 - to)p imt (x u x 2 ). (41) 

In the general case of correlated initial conditions, 
two additional terms containing contributions to 
the self-energy singular at to would contribute to 
the full G< see Paper III and Refs. [4,62]. 

Eq. (40) with all integrations explicitly shown 
reads 

G<(1,2) = 

J d3 J ^4G R (l,3)p< lt (3,4)G A (4,2) 

to to 

t\ t 2 

+ J d3 J ^/4G R (l,3)i: < (3,4)G A (4,2). 
to to (42) 

The first line involves, in fact, no time integrations; 
by (41), its structure is 

iG R (t u to)p imt G A (to,t 2 ). 



This term describes coherent evolution and gra- 
dual decay of the initial condition; for free 
particles, this would be the complete Gq . The 
other term, associated with the steady stage of the 
transport processes gradually takes over. In fact, in 
the long time limit to — ► — oo (...kinetic regime), the 
first term of (40) will die out and only the second 
term of the Dyson-Keldysh Equation will survive. 

While this integral form attributed to Keldysh is 
very illuminative concerning the causal structure 
of the theory and will be handy in formal 
manipulations, it is prohibitively complex for 
numerical work. A practicable form of a closed 
transport theory at the level of double-time non- 
equilibrium correlation functions is based on a 
(numerical) solution of the pair of coupled integro- 
differential equations, say (36) and (37) with the 
initial condition (39). The other propagator G A is 
found by conjugating G R . This is an alternative to 
the basic KB theory working in the same way with 
equations for G > , G < . 

7.5. Precursor equation 

To continue from the NGF to transport 
equations for the quantum distribution function, 
we use, instead of (37) or (38), their difference, 

{G R } _1 G < - G < {G A }- 1 = - G R I < + G A , 

(43) 

while retaining the Dyson equation for the 
propagators. The new equation is commonly 
called the precursor transport equation, because it 
already has a structure closely related to transport 
equations. Using Eq. (33), it can be rewritten with 
an unrenormalized drift on the l.h.s.: 

Gq 1 G < - G < Gq 1 = -G R I < +I ,< G A 

+Z R G < - G < Z A , (44) 

Employing (6), (34), we find the l.h.s. of the last 
equation at equal times: 

l.h.s. of (44) ^ ^ + i[H , p]_ . (45) 

This part of the equation has, thus, the form of the 
GME. The "interaction terms" on the r.h.s., how- 
ever, are expressed by means of the double-time 
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quantities. Our aim will be a reduction of these 
terms to expressions specified by a single time 
distribution function, p in the case of the GME. 

To obtain a kinetic equation involving the 
quasi-particle distribution /, we will have to 
proceed in a modified way. In particular, the 
drift term on the l.h.s. will be partly renormal- 
ized, that is, it will be somewhere between the 
fully renormalized drift in (43) and the bare drift 
in (56). 

Note on Dyson equations for g < 

There are reasons to favor one or the other 
convention as regards the correlation func- 
tions. Unfortunately, here we have to quote both: 
for historical reasons, use of g K , g > became 
common in papers concerning the quantum kinetic 
equations and the KBA. On the other hand, the 
GME, as a systematic tool preserving the innate 
structure of the NGF formalism, is more closely 
connected to the G < , G A choice. In agreement 
with this, we will use different conventions in 
Papers II and III. 

With the relations (compare Eqs. (5), (7) with 
Eq. (23) as discussed already in Section 7.1), 

G < (l,2) = i^ < (l,2), (46) 
G>(1,2) = -i^(l,2), (47) 

we have, for example, 

G R (1,2) = -i0(t 1 -t 2 )(g>(l 9 2) + g < (l 9 2)) 9 (48) 
G A (1,2) = i6(t 2 - t x ){g > {\ 9 2) + g < {\ 9 2)). (49) 

We recall that in the Langreth-Wilkins formal- 
ism, all equations are necessarily linear in the 
"less" quantities, while there is no such limitation 
on the propagator-related quantities. As a result, 
all previous results for G < , Z ,< are easily trans- 
formed to a form appropriate for g < . The 
correlation functions were related by (46), (47). 
Extending the same definition also to the self- 
energies and to the initial condition, that is, 



writing 

Z < (l,2) = i(7 < (l,2), (50) 

Z > (1,2)= -w>(l,2), (51) 

p<(l,2) = ip<(l,2), (52) 
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we obtain equations equivalent to (40) and (37), 
(38): 

g< = G R p< lt G A + G R c<G A , (53) 
and 

{GW = ( t < G a , (54) 
g < {G A }~ 1 = G R a < . (55) 
Finally, for the precursor equation we obtain 

-i(G« l g< -g < G^) =i(G R rr< -a<G A ) 

-i(l R g< -g<Z A ). (56) 



8. The reconstruction procedure 

The general plan of generating quantum trans- 
port equation within the NGF scheme is to take 
the precursor transport equation in one of its 
forms (43), (44) and (56), and to express the 
precursors to scattering (interaction) integrals in 
terms of one or the other distribution function. 
This will transform the equation involving two- 
point GF to a transport equation having the 
single-point distribution function as the unknown. 
For the purpose of this transformation, one of the 
approximate "Ansatzes" will be used, as outlined 
in Section 1. We wish to outline this approach in 
the present paragraph concentrating on the 
associated general questions. One basic question 
is: what is the real justification for such a round- 
about way to the transport equations? In other 
words, what is gained by working this way, 
compared to other approaches? 

Viewed from another angle, we may detach the 
construction of the transport equations from the 
fundamental issue that we already summarized as 
the hypothetical Reconstruction Theorem. To 
recapitulate the situation, it appears natural that 
virtually all of the relevant information about a 
non-equilibrium many-body system can be un- 
folded from its reduced characteristics, namely a 
pair of two-point quantities, say G < and G > . The 
Bogolyubov principle (8), just as the Ansatzes hint 
at the possibility that, actually, it may be enough 
to know or control even less, just a function of a 
single time variable, the one-particle density 
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matrix, or even a one-point quantity, the local 
particle density. Thus, the reconstruction problem 
of Section 3 can be made more specific: under 
which conditions the full description of the many- 
body interacting system can be built up from the 
knowledge of its single- time single-particle char- 
acteristics! 

It should be cautioned already here, however, 
that any Ansatz scheme in practical use, be it 
Kadanoff-Baym Ansatz, Generalized Kadanoff- 
Baym Ansatz, or any of their variants, plays the 
role of a physically well-justified, but still approx- 
imate truncation of the NGF. This would make 
the reconstruction procedures approximate al- 
ready by their nature. Another problem occurring 
in the Ansatz schemes is that the knowledge of the 
distribution function alone is not sufficient: the 
additional quantities entering the equations are 
propagators, and it has to be clarified, whether 
they are generated from the known distribution 
function in a separate step, or whether they are 
entering the scheme as an independent input, 
difficult to obtain beyond an approximation. Last 
but not least, we again have to caution that it is 
rather difficult to extend an Ansatz procedure so 
as to incorporate correlated finite-time initial 
conditions. These problems of principle have, of 
course, only a limited significance for practical use 
of the Ansatz technique, if it turns out to be 
convenient and gives reliable results. 

Let us have a look at the status of some 
alternative approaches vis-a-vis the reconstruction 
problem. Only then, will we be able to formulate 
the reconstruction hypothesis properly. 

First, the BBGKY method of an infinite chain 
of equations for reduced density matrices of all 
orders is inherently associated with a decoupling 
(truncation) in the chain [12,42-44,9]. This cer- 
tainly leads to a reduced description, but the whole 
approach is not directly related to the reconstruc- 
tion problem in the narrow sense followed here. 

Second, the projection methods in the Liouville 
space are a more likely candidate [64-70]. A 
partitioned expression for the projected (relevant) 
density matrix is closed and exact, in principle. It is 
open to generating ever-improving approxima- 
tions in a systematic way. The difficulty with this 
approach is its complexity, which exceeds that of 



the NGF method. There are some very important 
developments aiming at a combined technique 
permitting to take advantage from both sides, and, 
in particular, to extend the GF methods directly to 
systems with open boundary conditions at finite 
times [71-73]. 

There is a third, truly likely contender, the Time 
Dependent Density Functional Theory (TDDFT) 
[74-76]. In fact, some of the general results 
obtained in the context of this approach are about 
closest to the fundamental formulation of the 
Reconstruction Problem. At the same time, 
TDDFT draws from an independent source of 
physical inspiration. In view of this, we devote a 
paragraph to the TDDFT before discussing our 
main subject, the Ansatz approach. 

8.1. Functional substitution: inspiration by TDDFT 

A well-known case of functional substitution, 
seminal for the present topic, is the Hohenberg- 
Kohn Theorem concerning equilibrium many- 
body systems [76-78]. Their equilibrium state is 
specified by the time independent external field 
U(r) acting on the system. In particular, the local 
particle density n(r) depends on U(r). The two 
functions are dual: they enter the total energy 
<^ Hy of the system as the bilinear functional 
(U) = J d 3 rn(r)U(r). This suggests a functional 
substitution U n permitting to formulate the 
whole many-body problem as depending on the 
new functional variable. The theorem claims that 
such a substitution is possible — under conditions 
we do not care to elaborate. The new energy 
functional is obtained by the Legendre transfor- 
mation as E(n) = <C Hy — (U), which has been 
given on a firm mathematical basis [78,79]. 

As has been pointed out already in [17], " ... the 
feasibility of an exact building up of the full 
correlation functions from their time-diagonal 
sections can be motivated on physical grounds 
resembling somewhat the reasoning leading to the 
Hohenberg-Kohn theory." It became clear over 
the time that the non-equilibrium analogue of the 
Hohenberg-Kohn Theorem will be far more 
complicated: non-stationary dynamics of a many- 
body system is specified by two groups of physical 
constituents, the initial conditions and the driving 
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forces. Let us consider the simplest, and simulta- 
neously the particularly important case of the 
switch-on states. The first group, the initial 
condition, then concerns the system in equilibrium 
and includes the particle types and numbers, the- 
many-body interactions and the boundary condi- 
tions given by the thermodynamic environment. 
The second group, barring non-mechanical per- 
turbations, consists of only one, but crucial 
member; the external field U(x, t) acting from the 
initial time to on. This external field uniquely 
specifies each non-stationary mechanical process 
in the system. It enters the NGF formulation 
through the free Green's function G , Eq. (34). All 
quantities including the full Green's functions are 
functionals of the external field. The function dual 
to the space- and time local field is now given 
by the bilinear functional Jdln(l)U(l) as the 
particle density n(l) = p(l, 1). It is natural, but 
premature, to expect that the functional substitu- 
tion U —> n will also work in the non-equilibrium 
case, although this expectation has been a tacit 
assumption of the TDDFT for a long time. 

A firm basis for the subsequent development of 
the TDDFT is given by the Runge-Gross Theo- 
rem [80,74-76]. It plays the exactly role of the 
Kohn-Hohenberg Theorem, namely it is an 
existence theorem. It states that, for any initial 
condition specified by a pure many-body state fo 
at the initial time to, different external fields U(l) 
induce different local particle densities ^(1), so that 
the substitution U n is justified. (The fields U 
must be Taylor expandable in time at each time 
point). Here, we do not describe the specific 
problems of the TDDFT, connected to the time- 
dependent Kohn-Sham formalism of equivalent 
ghost systems of non-interacting particles. The 
Runge-Gross Theorem has a much wider applic- 
ability, however. 

One nice point is that the Theorem provides a 
proof of the Bogolyubov conjecture quoted in 
Section 3, at least for a class of typical physical 
situations, which is not overly specialized. In fact, 
the Theorem implies even more, because the 
reduction to the one-particle distribution applies 
right from to, while the Bogolyubov postulate 
allows for a formation time x c which must elapse 
before the reduction takes place. This notion is 



physically appealing, and it would be interesting to 
study the early Runge-Gross evolution of the 
system in more detail. 

It should be stressed that an entirely arbitrary 
initial condition cannot be uniquely characterized 
by its one-particle reduced characteristics. This is 
reflected in the formulation of the Theorem, of 
course. We may note that the situation is similar to 
the one we have already encountered in connection 
with the correlated initial condition for the NGF. 
There, we have singled out three special classes of 
simple initial conditions: the uncorrelated initial 
condition at to -> — oo (the Bogolyubov decay of 
initial correlations), the finite to initial condition 
uncorrelated by assumption (Keldysh IC) and the 
switch-on states. The last case is particularly 
appealing: the initial state is generated by a 
perturbation of a "dark" (i.e., no field) state in a 
distant past. If this dark state is a state of 
equilibrium, the Hohenberg-Kohn Theorem ap- 
plies; then, the state evolves under the conditions 
of the Runge-Gross Theorem. Consequently, here 
we encounter an initial condition fully reducible to 
the history of n(r, t) for t preceding to. 

The last point of this section focuses on the 
actual implementation of the functional substitu- 
tion. The original idea to base it on the real time 
Legendre transform jd\n(X)U(X) encountered 
inconsistencies, as mentioned above. As could 
have been expected, an adequate description of the 
non-equilibrium dynamic states calls for the use of 
the Schwinger-Keldysh trajectory, so that the 
proper modification of the functional becomes 

J dl/!(l)tf(l) = <y | J dlw(l)l7(l)|!Po> 

c c 

es i J dihuy. 

c 

Recognition of this integration path was one of the 
key steps in developing the contemporary TDDFT 
and linking it with the NGF language and 
technique [75,81]. For all details on this extensive 
and developing field, the reader is referred to 
reviews and to original papers. It is rewarding to 
observe that by this development, two major 
branches of the Schwinger school, the DFT as 



168 



V. Spicka et al / Physica E 29 (2005) 154-174 



represented by W. Kohn, and the NGF with the 
Schwinger-Martin-Kadanoff & Baym/Horing 
pedigree [21,82,81], converge together in one major 
contemporary direction of many-particle/con- 
densed matter theory. 

8.2. Intuitive way to reconstruction: 
Kadanoff- Baym Ansatz 

To physically motivate the general formulation 
of the reconstruction problem in the context of the 
NGF based transport theory, we will introduce the 
historically first solution of the reconstruction 
problem and stress some of its important features. 
This was proposed by Kadanoff and Baym in 1962 
[3], when they aimed at the construction of a 
quantum Boltzmann equation from the KB 
equations for g < and g > . Their approach thus 
did not strive for generality. They rather made two 
appropriate assumptions: (i) The disturbance U(l) 
had to vary slowly both in space and time, (ii) The 
excited states engaged in the transport process 
were to be stable quasi-particles. 

Here, we will not follow the original work, but 
rather show, without detailed references, how this 
line of thinking developed over time. The two 
assumptions are formalized in two steps. First, the 
transition to the Wigner representation for all four 
space-time variables is made. For any two-point 
function, F(l, 2), we make the transformations 



[x u til [x 2 ,t 2 ] 



F(xi;tux 2 ;t 2 ) 



[r = 
[x 



1 

2 

: X\ 



(xi +x 2 ),t = Uh + f 2 )], 



■x 2 ,t = h 



tl] 



F.T. 



F(r + -x;t- 
f(k, co, r, t) 



-t,r — -x;t ■ 



in 

(57) 



This leads to the notion of a local electron 
structure, slowly dependent on the changing 
average position r and time t. In the second step, 
Eq. (55) for g K is written in the Wigner (mixed) 
representation and all damping terms are left out, 
accounting for the assumption of stable quasi- 
particles: 

[co-e k - U eff (r, t) - a(k, co, r, t)]g K (k, co, r, t) = 0. 

(58) 



We use the shorthand a — Reo" R , see Eq. (11.14) in 
Paper II. The general solution of this equation is 

g < (k,co,r, t) = (j)(k,to,r, t) 



xd (co — Ck — U e ff(r, i) — a{k, co, r, t)) , 



(59) 



where the second ^-factor allows in the present 
case to define the local quasi-particle dispersion 
law by 

s(k, r, t) = c k + U eff (r, t) + o(k, co, r, t)\ m=<kAty 

(60) 

while the first factor is an arbitrary function. 
Finally, solution (59) assumes the form, which will 
presently be termed the Kadanoff- Baym Ansatz: 

g < {k, co, r,t) —f (k, r, t) x S(co — s(k, r, t)), 



f(k,r, t) = 



1 



1 9_£ 

1 6co 



cj)(k, co, r, t) 



(61) 



co=£(k,r,t) 



Here, 3(co — s(k, r, t)) obviously plays the role of 
the quasi-particle spectral function Aq(co,k, r, t). It 
is then natural, and later will be verified, that the 
function / should be identified with the quasi- 
particle distribution. In Paper II, this identification 
will be made by employing (61) in the precursor 
equation (56) and demonstrating that / obeys the 
quantum kinetic equation. 

Here, we link (61) with the Fluctuation Dissipa- 
tion Theorem rewritten for g K in the form 



g (k,co,r,t) = g (k,co) 



= ./fdM x A(k,co). 



(62) 



This identity is strictly an equilibrium feature and 
reduces the number of independent Green's func- 
tions to one: /fdO^) ^ s tne Fermi-Dirac factor, a 
univeral function of energy. If A is replaced by ^q, 
we obtain the equilibrium limit of (61) with 



f(k,r,t)=f(k)=f FD (e(k)) 



(63) 



This is precisely the equilibrium quasi-particle 
distribution. We see that the Kadanoff-Baym 
Ansatz can be viewed as a non-equilibrium exten- 
sion of the Fluctuation Dissipation Theorem to 
locally equilibrated states, suitable not very far 
from equilibrium. In qualitative terms, beyond 
equilibrium, the two functions g < and A are not 
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completely reduced one to another: additional 
information, the quasi-particle distribution /, is 
required. Still, g < plus A can be reconstructed from 
a substantially reduced information contained in / 
plus A. 

To gain a better feeling of the cavalierly 
interchanges A < — > Aq, let us introduce the 
renormalization constant z and rewrite the exact 
equilibrium particle spectral density as 



A(k, co) = z(k)d(co - s(k)) + Satellites 



z(k) 



1 



1 



8co 



(64) 



co—s(k) 



Let us try to leave out the satellite contribution. 
This permits to rewrite (62) in two "improved" 
KB A- like ways: 



= / x z x 3(co — £)==< 



QP 



particle ^4 



x z • <5(co — £) 
f - z x 5(co — s) . 

QP 



(65) 

Here, the particle correlation function g < fac- 
torizes either to the QP distribution times the true 
spectral function (upper line), or to a product of 
the particle distribution and the QP spectral 
function (lower line). In either form, a most 
important aspect of the true Kadanoff-Baym 
Ansatz is lost, namely an automatic particle 
number conservation. This is achieved, in the 
KBA, by lumping together the satellite fraction 
1 — z of the spectral weight with the z weight of the 
quasi-particle resonance. Clearly, a better treat- 
ment of the satellites is desirable. This leads to the 
concept of the extended quasi-particles explained in 
Paper II. 

Another important comparison between the 
Fluctuation Dissipation Theorem (63) and its 
KBA extension (61) is best made in the original 
space-time representation (only the time argu- 
ments are shown for clarity): 

G < (t u h) =f( t ^^)[G A (t u t 2 ) - G R (t u h)l (66) 

where the spectral function is written as 
A = i(G A — G R ), see (30). In correspondence with 



the quasi-classical nature of the KBA, the quasi- 
particle distribution is taken at the mid-point time 
t = \{t\ + t 2 ). This feature violates the overall 
causal structure of the NGF theory manifested in 
Eqs. (40) and (42). If the Ansatz is to consistently 
extend the distribution function to G < , the 
distribution should appear as an initial condition 
of its kind at the earliest time and then be 
propagated to later times. In Section 4 of Paper 
II we will discuss how the acausal features of the 
KBA are closely related to the description of the 
particle collisions in the quantum BE. The KBA 
time structure (66) becomes very disturbing 
beyond the quasi-classical limit, and this was one 
of the reasons for introducing the so-called 
GKBA, which was constructed to be causal from 
the outset. This Ansatz is the subject matter of 
Paper III. The basic ideas leading to the GKBA 
and from there to the "Reconstruction Theorem" 
will be outlined, in the next subsection. 



8.3. From GKBA to the Reconstruction Theorem 

The GKBA has been introduced as an alter- 
native to the original KBA [17-20,8,9]. It can be 
written as 

G < (t u t 2 ) = -G R (t u t 2 )p(t 2 ) + p(ti)G A (t u t 2 X 

(67) 

In comparison with the KBA, as given by (66), it 
differs in terms of four constitutive properties: 

Causal structure: The two propagators are now 
separated and each is multiplied by its 
own distribution function corresponding 
to the earlier time argument, depending 
on the retarded or advanced propagation 

Equal time limit: The distribution function is the 
true one-particle density matrix, not the 
quasi-particle distribution. Thus, at equal 
times, the GKBA for G < obeys identity 
(6). This implies an automatic particle 
number conservation. 

General representation: The order of factors cor- 
responds to the products interpreted 
as a matrix multiplication. Therefore, a 
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diagonal representation for Green's func- 
tions is not assumed. 
Arbitrary non- equilibrium: No assumptions about 
slowly varying disturbances, small devia- 
tions from equilibrium, or quasi-particle 
approximation are built in, the Ansatz is 
formally quite general. 

The GKBA particle correlation function perfectly 
fits into the formal NGF scheme. It was con- 
structed in this manner in order to overcome the 
limitations of the KBA. When inserted into the 
machinery of deriving quantum transport equa- 
tions from the NGF, it leads, as described already 
in the Introduction, to the GME, rather than to 
the BE like quantum kinetic equations. For these 
purposes, it has been amply used. All this is the 
subject matter of Paper III. 

The GKBA also has an important meaning, 
independent of and going beyond the restricted 
context of transport equations. Namely, it offers 
an alternative path to the Reconstruction Pro- 
blem, as will be outlined here. 

Comparing (67) with (62), we may state that the 
GKBA is the utmost extension of the Fluctuation 
Dissipation formula to non-equilibrium. The related 
interpretation resembles that of the KBA: starting 
from the reduced information in p and the 
propagators, the GKBA permits to reconstruct 
the whole particle correlation function as their 
product. 

There is a basic difference, however, which 
distinguishes the GKBA from the KBA. Namely, 
while both Ansatzes are but approximate con- 
structions, expression (67) can, in addition, serve 
as a starting approximation for a process leading 
to an exact reconstruction of G < . 

This assertion is valid in this simple form, 
and can be demonstrated, for the case of 
uncorrected initial conditions. From Eq. (40), 
there follows a couple of appropriate Dyson 
type equations for G < , reported in [17] and 
discovered independently by Vinogradov [18]; see 
also [8,9]. 

To illustrate their derivation and interpretation, 
we start by observing that the GKBA expression 
(67) for G < consists of two distinct parts, one for 
t\ > t 2 , the other one for t 2 > t\ and The exact G < is 



split in the corresponding way: 

G < (t x , h) = G< (h , ti) - Gt (t x , t 2 ), (68) 
G<(t u t 2 ) = 6(t l -t 2 )G < (t u t 2 l (69) 
GZ(t u t 2 ) = - 0(t 2 - t^ituh) = [G<(t 2 ,h)f . 

(70) 

It is convenient to calculate 

{G R r l G<\ tut2 

= (Go 1 - Z R )G^\ tut2 
= S(h - h)p(t 2 ) 

+ 9(t l -t 2 )({G R r l G < 

+ Z*G< -Z R G<)\ tut2 
= 8(ti - t 2 )p(t 2 ) 

+ e(t l -t 2 )(z < G A + z; R G<)\ tut2 . (7i) 

where Eqs. (33), (34) and (37) were employed. 
Multiplication by G R from the left yields the first 
equation for G^, that is for G < valid in the 
t\>t 2 > to wedge: 

G < (t u t 2 )= -G R (t u t 2 )p(t 2 ) 

+ f 1 dl [ 2 dtG R (t l J)I < (l t)G\l t 2 ) 

Jt 2 J t 

+ [ 1 df / 2 dtG\t x J)Z\l t)G < {t, t 2 ). 

Jt 2 J t 

(72) 

In the complementary time region t <t\<t 2 , the 
equation for G < = G^ is quite similar. It can be 
obtained as a conjugate to (72), as it holds that 
G < {tiJ 2 ) = -[G < fe,^i)] t . It starts as 

G < (t u t 2 ) = p(t l )G A (t u t 2 ) 

+ • • • . (73) 

The absolute terms of equation pair (72) and (73) 
together form the GKBA expression (67). The two 
integral terms contain the unknown G < : the 
second one explicitly, and the first one through 
the functional dependence I K = I < [G < ]. 

We may conclude that by solving Eqs. (72) and 
(73), the double-time correlation function 
G < (t\,t 2 ^t\) is generated from the knowledge of 
its time-diagonal part, the one-particle single-time 
density p(t) oc G < (t,t), and of the propagators 
also taken as known. This result in itself has the 
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importance of the exact non-equilibrium counter- 
part to the Fluctuation Dissipation Theorem. 
Furthermore, it represents the first, and the most 
important step to the NGF Reconstruction 
Theorem. 

Eqs. (72) and (73) can also be used as they stand 
in two alternative ways. On the one hand, they are 
perfectly suited for a direct numerical solution 
aiming at G < connected with a particular physical 
situation. So far, no computation of G < along 
these specific lines has been reported. On the other 
hand, the reconstruction equations (72) and (73) 
can serve the important purpose of generating 
gradually improving Ansatzes by formal iteration. 
The use of such an improved renormalized GKBA 
Ansatz would be to introduce it into the interac- 
tion terms of the precursor transport Eq.(56) in 
order to convert it to a closed GME for p. This is a 
well-defined route to a family of presumably more 
and more precise GME. 

8. 4. The NGF Reconstruction Theorem 

To clarify the meaning of the Theorem, first we 
have to see, what is the necessary input for 
obtaining G < . (i) At equal times, Eqs. (72) and 
(73) turn to tautology. The "time diagonal" p is 
one of the input quantities and has to be known 
beforehand. (This trivially includes the un-corre- 
lated initial condition (39).) (ii) We assumed 
throughout that the propagators and their self- 
energies are somehow known. This is not trivial 
beyond equilibrium and, in fact, the propagators 
have to be found in another process, too. 

The NGF Reconstruction Theorem is obtained 
by choosing an appropriate order of steps in 
solving the full NGF problem. At the hub of the 
whole process stands the single-particle distribution 
p. It is obtained by solving the exact precursor of 
the GME, the so-called GKB Equation, which is 
obtained by joining Eqs. (44) and (45) 

^ + i[H ,p]_= -i[-G R I< +I<G A 
ot 

+ S R G< -G<Z\ =t=t2 . (74) 

To convert the GKBE to a closed transport 
equation for p, we have to specify the physical 



ill 

approximation for the self-energies in the form 
I = E[G], that is, more explicitly, 

(75) 

Next, the pair of reconstruction equations (72), (73) 
is used to express G < in the functional form 

G K = G < [p,G R ' A ] . (76) 

Finally, the Dyson equations for propagators (see 
(36)) 

G R = G R + G R Z R G R , 

g a = g a + g a z a g a . (77) 

link the problem with the non-interacting particles 
and incorporate the effective external field as a 
boundary condition. 

The full NGF set G R , G<, G A is obtained by a 
simultaneous solution of Eqs. (74)-(77). 

The NGF Reconstruction Theorem 

The contents of the Reconstruction Theorem 
appears to be such that the determination of the 
full NGF set G R , G < , G A can be structured to a 
dual process. In one stream, the single-particle 
distribution p is obtained from a closed GME type 
transport equation. In the other stream, the NGF 
is reconstructed from known p. The two streams 
are mutually coupled: 



p 


— > 


G 


GME 




reconstr. 




<- 





(78) 



This prompts a simultaneous treatment of the 
whole task, which is only possible by a solution in 
successive steps. 

These steps may progress incrementally in time. 
Starting from the initial p(to) = p init , see Eq. (39), 
all quantities are obtained step by step rimming 
the already covered time region [to, t] (g) [to, t] by a 
strip having a width of the time increment At. It is 
easy to verify that all equations involved have the 
causal structure appropriate for this protocol. The 
experience with similar tasks suggests that this will 
be a practicable technique [62]. 
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An alternative method would be an iteration of 
the equations with respect to some expansion 
parameters. While less promising for computa- 
tions, the iteration scheme has a more basic 
context. Namely, it is important to select and 
interpret the iteration procedure. There are two 
standpoints here, and their relationship is not 
entirely clear. 

/. Interaction strength as a small parameter 

Clearly, for non-interacting particles, the self- 
energies and Eqs. (72) and (73) demonstrate that 
the GKBA is exact in this limiting case. Thus, the 
iteration of the equations might be interpreted as a 
perturbation expansion depending on the particle 
interaction strength. This point of view has been 
expressed several times, and it is not unjustified. If 
the whole reconstruction reduced just to this, it 
would be brought rather close to the direct 
methods of deriving the GME, which are often 
based on some notion of a weak influence of the 
particle interactions. 

//. Collision duration time as a small parameter 

When comparing the original Eq. (42) for G < 
with (72) and (73), an important characteristic fea- 
ture of the latter is conspicuous: the off-diagonal 
integration regions of the double integrals: 

Both integration variables converge in the self- 
energies I < (t,i), I R (t,i). Only the off-diagonal 
blocks of the self-energies joining the past and the 
future are contributing. The small parameter 
determining the importance of the corrections to 
the GKBA in, say, Eq. (72) is, loosely speaking, 
the width of the strip around the time diagonal, in 
which the values of self-energies are significant. 
The respective characteristic times, tq for I R and 
t c for I ,< , may be different, but both should be 
"small". On the physical side, we see that, rather 
than the strength of the interactions, it is their 
inner dynamics which is decisive for the iteration 
procedure. This dynamics is reflected in the 
memory for build-up processes induced in the 
system during a non-equilibrium evolution. In 



fact, the two characteristic times are best defined in 
an operational way from Eqs. (72) and (73). 
Whether they are small or not, that depends on 
the process in question and the related character- 
istic times. Thus, there is no absolute criterion for 
validity of the GKBA or its iterates. 

It should be pointed out that with this second 
interpretation of the Ansatz scheme, the whole 
path to the quantum transport equations outlined 
in this section has practically the same physical 
background as the direct methods based on the 
superprojection techniques. There may be a 
practical difference, however, in the relative ease 
of actual construction of the quantum transport 
equation in the NGF method, where the elements 
used are all single-particle, although double-time, 
characteristics and the whole scheme is consistently 
renormalized at each level of approximation. 

Another comparison relates to the TDDFT. The 
NGF Reconstruction Theorem concerns the stage 
of reconstruction up from the one-particle density 
matrix. The Runge Gross Theorem is more specific 
in showing that the necessary reduced quantity is, 
in fact, the local particle density, that means only 
the space diagonal of the density matrix. This is 
correct for a local external field; in our treatment, 
no specific property of the external potential was 
employed, in fact. For the more general situation 
of a possible non-local field, the reduction to the 
full p is to be expected. The other difference 
concerns the proof of the Theorem. In the NGF 
case, it is constructive, so that the proof also 
contains the algorithmic basis for an actual 
reconstruction. Otherwise, the relationship of both 
approaches remains close, and hybrids, or off- 
springs of both may be promising in the future. 

9. Conclusions 

This paper served several purposes: 

(1) to summarize the present views on construc- 
tion of the quantum transport equations based 
on the NGF approach; 

(2) to characterize the position of the NGF among 
other theoretical techniques for non-equili- 
brium systems; 
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(3) to analyze the physical and formal context of 
the so-called Ansatz decouplings; 

(4) to explore the status of the so-called Recon- 
struction Theorems, which try to show that an 
exact description of a non-equilibrium system 
in terms of one-particle quantities is possible. 
Such ideas are at the core of the density 
functional theory. We presented a comparison 
between the TDDFT and the properly renor- 
malized NGF formalism. There is a close 
relationship between both formalisms and they 
may tend to merge in the future. 

The general outcome of our analysis is favor- 
able. The discussed approaches and their results 
have been rather successful and form a basis for 
future work aiming at a description of complex 
processes far from equilibrium, typically involving 
coherent short time dynamics in response to fast 
disturbing fields, which is followed by early 
relaxation and equilibration processes evolving 
on the intermediate time scale. 

The two attached Papers II and III are devoted 
to a more detailed treatment of the different time 
scales for quantum transport processes and to a 
lesser extent to the bridging periods between them. 
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